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Abstract 

We study the dynamical fluctuations of various particle yield ratios at different 
incident energies. Assuming that the particle production yields in the hydronic final 
state are due to equilibrium chemical processes (7 = 1), the experimental results 
available so far are compared with the hadron resonance gas model (HRG) taking 
into account the limited momentum acceptance in heavy-ion collisions experiments. 
Degenerated light and conserved strange quarks are presumed at all incident energies. 
At the SPS energies, the HRG with 7 = f provides a good description for the measured 
dynamical fluctuations in + /ir~)/(7r+ + tt~). To reproduce the RHIC results, 7 
should be larger than one. We also studied the dynamical fluctuations of (p+p)/(7r+ + 
7r~). ft is obvious that the energy-dependence of these dynamical fluctuations is non- 
monotonic. 

1 Introduction 

Understanding of the dynamical properties of hot and dense matter is a key question 
in the heavy-ion collisions experiments. The phase structure and event-by-event fluctua- 
tions [H El [3] have been suggested to provide comprehensive characteristics of the particle 
production yields. They are essential observations to examine the hypothesis about the 
equilibrium of the chemical processes in the hadron final state j4] . The event-by-event fluc- 
tuations of certain particle yields have been studied at SPS and RHIC energies |5l El E]- 
Therefore, it is natural to study the energy-dependence of the particle yield ratios and the 
event-by-event fluctuations using the hadron resonance gas model (HRG), as it provided 
a good description for the thermodynamical evolution of the hadronic system below the 
critical temperature [51 El [TU] and has been used to characterize the conditions deriving 
the chemical freeze-out |111 [T2] . 

The hypothetical chiral symmetry breaking restoration and deconfined phase transition 
to quark-gluon plasma (QGP) are to be characterized by remarkable fluctuations in the 
particle production yields [fl El [3], which are accompanied by dynamical and volume 
fluctuations as well. The latter can simply be eliminated, when taking into consideration 
the dimensionless particle yield ratios [3]. 

In this letter, we try to answer the questions, whether the strangeness quarks Qs should 
enhance the dynamical fluctuations and whether the critical endpoint could be localized 
by means of event-by-event fluctuations in the hydronic final state. We make predictions 
for the dynamical and statistical fluctuations of the ratios of different particle yields in 
dependence on incident energy. Apparently, the dynamical fluctuations strongly depend 
on sort of the particle yields. In some particle yield ratios, the dynamical fluctuations are 
smaller than the statistical ones. In others, the dynamical fluctuations are slightly greater 
than the statistical ones. The energy dependence is non-monotonic. 



*drtawfik@mti.edu.eg 



All these predictions are phenomenologically of great interest and hope to encourage 
experimental attempts to measure their event-by-event fluctuations in a wide range of 
incident energy. 



2 Model 

The hadron resonances treated as a free gas [U EJ [lOl [El [Hj are conjectured to add to the 
thermodynamic pressure in the hadronic phase. This statement is valid for free as well 
as strong interactions between the hadron resonances themselves. li has been shown that 
the thermodynamics of strongly interacting system can be approximated to an ideal gas 
composed of hadron resonances [HI [15] 

At finite temperature T, strangeness and baryo-chemical potential /x^, the pressure 
of one sort of hadron resonance reads 



1 ± 7 A^Ase T 



e(k) 



(1) 



where e{k) = {k"^ + vn?'Y^'^ is single-particle energy and it stands for bosons and fermions, 
respectively, g in the front of the integration is the spin-isospin degeneracy factor and 
7 = 7^7™ stand for the quark phase space occupancy parameters, where n and m being 
number of light and strange quarks, respectively. A = exp(^/T) is the fugacity, where n 
is the chemical potential multiplied by corresponding charge. Summing over all hadron 
resonances results in the final thermodynamic pressure in the hadronic phase, as no phase 
transition is conjectured in HRG. 

The quark chemistry is given by relating the hadronic chemical potentials to each 
of the quark constituents; = and fis = fJ^q — fJ-s, where q and s being the light 
and strange quark quantum number, respectively. The baryo-chemical potential for the 
light quarks is averaged as fig = -|- fid)/'^ and the strangeness chemical potential 
fis is calculated as a function of T and fiB under the assumption of strange quarks 
conservation 1131. 



In grand canonical ensemble, the particle density is no longer constant. 

i 

{{ANf) = yA^fdke ^^'^ „^ (3) 

To count the number density N and fluctuations (AA^)^ in the hadronic final state, the 
chemical freeze-out processes have to be taken into account, i.e. the hadron resonances 
should finally decay to stable particles or resonances. 

(ATf-^) = {Nf^^^')+Y,bj->rm, (4) 

((AA,^,)2) = h,^i{l-h,^,){N,) + h]^,{{AN,f) (5) 

where hj^i is the branching ratio for the decay of j-th to i-th particle. In this work, the 
chemical freeze-out is characterized by s/T^ [14J, where s is the entropy density. 
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The fluctuations of particle yield ratios of particle 1 and particle 2 read [3] 



In this expression, we include all possible fluctuations, i.e., dynamical and statistical 
as well. The third term of Eq. [U] counts for the fluctuations from hadron resonances 
that decay into particle 1 and particle 2, simultaneously. In such a mixing channel, all 
correlations including the quantum statistical ones are taken into account. Obviously, this 
decay channel results in strong correlated particles and correspondingly fluctuations. 

To extract the statistical fluctuations, we apply Poisson scaling to the mixed decay 
channels 0, 

Detector acceptance factor and resolutions are main sources for statistical fluctuations. 
Subtracting Eq. [7] from Eq. [6l we get the dynamical fluctuations of the particle yield ratio 
iVi/A^2. 



3 Results 



The experimentally measured dynamical fluctuations of particle yield ratios 
(K~^ + K~)/(7r~^ + 7r~) are systematically confronted with the theoretical predic- 
tions in Fig. [H An earlier attempt to compare with preliminary results has been reported 
in [21 [5]. It was found that the theoretical and experimental ratios of dynamical to 
statistical fluctuations are compatible with each other. Individual fluctuations themselves 
are not. 

Depending on 7, HRG is apparently able to predict various particle yield ratios at a 
wide range of incident energy. At SPS energies, HRG with 7 = 1 provides a good descrip- 
tion for the experimentally measured dynamical fluctuations [5l[7]. To reproduce the RHIC 
results, 7 should be larger than one. The dynamical fluctuations of {p + p)/(vr"'" + tt") are 
depicted in Fig. [2j Few comments are in order at this moment. 

• The dependence on ^/s is non-monotonic. The fluctuations can be suppressed and/or 
enhanced at different 

• Strangeness fluctuations are positive and enhanced with ^/s. There are remarkable 
minima at the top SPS energies 

• At high energies, the fluctuations smoothly increase with y/s 

In Fig.[Tl we compare measured dynamical fluctuations of {K~^ + K~)/(tt~^ + it~) yield 
ratios with HRG model. At 7 = 1.0, we find an excellent agreement at SPS energies. At 
RHIC energy, the measured fluctuations are above the theoretical ones. One has to allow 
7 to have values large than one, in order to reproduce the data. The explanation for this 
disagreement would be two-fold. First, the RHIC measurements are still preliminary [7]. 
The final measurements might modify the 7-value reported here. Second, we refer to our 

^Experimentally, there are various methods to eliminate the statistical fluctuations [7]. The frequently 
used one is the counting of particle yield ratios from mixing events. 
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Fig. 1: The dynamical fluctuations of (K^ + K^)/ (vr^ + vr") ratio as a function of square 
root of center of mass energy ^/s (curves) compared with experimental results (circles). 
SPS results [5] are drawn in open symbols while preliminary RHIC results [7| are given by 
solid circles. HRG with 7 = 1.0 can very well reproduce the SPS results. Larger values 
have to assign to 7 in order to reproduce the RHIC dynamical fluctuations (see the text). 

previous study of the particle yield ratios in heavy- ion collisions [16], where we concluded 
that the statistical models, like HRG, with 7 = 1 slightly overestimate the particle yield 
ratios at RHIC energy. Therefore we should assign to 7 values larger than one. What we 
observe here apparently supports such a conclusion that the quark phase space occupancy 
factor, characterized by 7, is obviously modified at RHIC energies. A connection between 
7 and hypothetical phase transition is discussed in [T6]. Therefore, the modification of 
7-value might be understood, if phase transition to QGP takes place. 

Additionally, we see in Eq. [8] that {Ni), where i equals 1 or 2, inversely proportional 
to fj, assuming that the contribution from the two-particle-channel is negligible {N1N2). 
Should HRG model overestimate the particle yield ratios at RHIC energy, as the case 
studied in [16j, straightforwardly, we expect that the dynamical fluctuations a should be 
underestimated at 7 = 1.0. 

So far we can conclude that the quark phase space occupancy factor obviously depends 
on the incident energy. Up to the top SPS energy, the quark phase space occupancy 
factor is most probably saturated (7 = 1.0). At RHIC energies, a phase transition to 
new state of matter has been reported [IT]. This can be seen in our analysis, if the 
quark phase space occupancy factor is modified. Allowing 7 to have values larger than 
one, HRG bests reproduce the particle yield ratios as reported in |TT1 [12l [16] and their 
dynamical fluctuations, the subject of this work. 

The statistical model analysis of {K'^ +K~) / {tt^ +tt~) dynamical fluctuations reported 
in [18] let to the conclusion that the non-equilibrium model, equivalent to 7 7^ 1, would 
provide an acceptable description at top SPS and RHIC energies. The definition of the 
chemical freeze-out and the mechanism assuring conserved strangeness in the hadronic 
final state would be responsible for this contradiction with HRG. 
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At low incident energy, the energy density or temperature and the degrees of freedom 
are not high enough to cause non-equihbrium phase transition to QGP [13^ I14j. As the 
hadronic matter goes into a new phase of nearly deconfined quarks and gluons, the quark 
phase space occupancy factor should correspondingly change, The fitted values of 7 used 
to draw the figure 1 in [18] are not explicitly given to discuss them here. Taking into 
account the momentum acceptance, the HRG results show that 7 = 1 best reproduces 
the dynamical fiuctuations at low SPS energies. 

As the HRG calculations are performed in grand canonical formulism, the so-called 
conservation laws |19j would be taken into account. With the conservation laws we mean 
- among others - the deference between the canonical (CE) and grand canonical ensembles 
(GCE). At high energies, equivalent to large volume V, - in principle - there is almost no 
difference |19j between CE and GCE. The conservation laws are discussed at vanishing |19j 
and finite chemical potential in the final particle multiplicity [20]. In GCE, the control 
variables are fi, V and T. The total number density is therefore allowed to fiuctuate. 
Therefore N is related to CE by the Legendre transformation. 

In HRG, we studied this problem from another point of view [13] and found that 
the ideal quantum Boltzmann gas formulism applied in HRG takes into account the in- 
teractions and the quantum statistics in the system. Furthermore, the ideal quantum 
Boltzmann gas overcomes the limitations in classical Boltzmann ideal gas, that the en- 
tropy for instance is only specified within an undetermined additive constant in the way 
that it determines such an additive constant in the high temperature limits of quantum 
Fermi- and Bose-gas. Quantitative estimation of the conservation laws at finite chemical 
potential, i.e. different incident energies, is included in the parameter q [19], where q is 
probability to detect n particle density out of total N particle density produced in the 
whole momentum space. The averaged particle density detected (accepted) read 



where (N) and {{AN)'') are given in Eq. [2]- Eq. [5l 

In Fig.O the dynamical fluctuations of non-strangeness (p + p)/(vr^ -|-vr~) yield ratios 
are depicted as a function of ^/s. At q = 1, HRG obviously overestimate the dynamical 
fiuctuations. It is worth noticing that crayn excursively increase at y/s ~ 15 GeV. The 
negative values are to be understood as dominant statistical fiuctuations, especially the 
ones not explicitly included in Eq. [7l The negative values might also refer to dominant 
fiuctuations in the proton-pion channel. 

Taking into account an acceptance factor q < 1, the experimental data can be repro- 
duced. As noticed in previous figure (Fig. [2]), the quark phase space occupancy factor 
is expected to be modified at incident energies higher than SPS. Taking this fact into 
account, the dynamical fluctuations of {p + p)/(vr+ -|- ■k~) explosively switch to positive 
values. If this expectation turns to be correct, we expect that the dynamical fluctuations 
of {p + p)/{tt^ + vr~) at RHIC energies might not be as smooth as the SPS ones. 

4 Discussions and conclusions 

At high temperatures, the hadronic matter is conjectured to go through chiral symmetry 
breaking restoration and deconfined phase transition(s) to partonic phase at almost same 
critical temperature [21]. Depending on order of phase transition, which in turn depends 
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Fig. 2: The dynamical fluctuations of non-strangeness yield ratio (p + p)/(vr''' + vr^). The 
experimental results (open circles) are taken from [5l [6]. The negative values point to 
dominant statistical fluctuations. The factor q relates canonical with grand canonical 
ensembles (see the text). 

on quark flavors and masses, fluctuations in particle yields are likely expected, especially 
if the transition causes out-of-thermal equilibrium. The three pions, vr*^'^, are the lightest 
Goldstone bosons resulting from chiral symmetry breaking. Therefore, strong fluctuations 
in the pion flelds are likely expected during the transition in chiral symmetry. 

The chiral transition associated with massless up and down quarks is of second order. 
If the quark masses are accounted with, the transition is either a slightly flrst order or 
just a cross over. The transition is explicitly a first order, if strange quark with very light 
mass is included. Also the deconfinement phase transition depends on the quark fiavors 
and their masses. For three quark flavors, the transition is likely of flrst order. In this 
case the transition occurs via nucleation of hadronic bubbles in the background of QGP. 
For two flavors, the transition is a smooth cross-over. 

This discussion might illustrate how strong are the pion dynamical fluctuations con- 
nected with the chiral phase transition at high temperature and out-of-thermal equilib- 
rium. These fluctuations are given the generic name of disoriented chiral condensates 
(DCC) [22]. 

The dynamical fluctuations associated with strong flrst order of phase transition are 
likely very large. The continues second order or cross over phase transition might wish out 
large part of dynamical fluctuations in the flnal state. On the other hand, the dynamical 
fluctuations are conjectured to slow down near the second order phase transition. This 
has been confirmed in classical systems, solid state physics. In quantum field theory, 
the long-wavelength (spinodal) modes will be quenched through the second order phase 
transition [23j. 

The Fluctuations of quark number have been studied in lattice QCD [24J. It has 
been found that the T-dependence of Uq fluctuations is dominated by the analytic part 
of the partition function. Across the critical temperature, there is a smooth increase in 
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the fluctuations. Almost same results are reported in this work. We therefore can use the 
fluctuations to characterize the phase transition. 

In our analysis, the critical temperature is not exactly specified, as HRG does imple- 
ment any phase transition. The most economic way is to study the dynamical fluctuations 
in the hadronic final state, such a way we can also compare with the experimental results. 
We picked up two particle yield ratios, for them we have extensive experimental estima- 
tions at a wide range of incident energies. We find that the dynamical fluctuations depend 
on the particle yields and incident energy. 

TabU] gives the relation between dynamical and statistical fluctuations Udyn/ (^stat at 
SPS and RHIC energies. For the strangeness particle yield ratios, the statistical fluctu- 
ations increase against the dynamical ones. For non-strangeness particle yield ratios, we 
get just the opposite. 



0'dyn/(^stat 


12.3 


17.3 


62.4 


200 


{K + K-)/{'K+ + TT-) 


0.0143 


0.0140 


0.0115 


0.0088 


(p + p)/(7r+ + 7r~) 


-0.0143 


0.0139 


0.0309 


0.0317 



Tab. 1: ct dyn / <^ stat of the particle yield ratios at 7 = 1 and q = 0.75 and various y/s. Using 
this quantity, we can estimate how large are the dynamical fluctuations compared to the 
statistical ones when the produced particles are chemically frozen. The validity of these 
predictions apparently depend on the conditions controlling the chemical equilibrium in 
the flnal state. 

The energy dependence of dynamical fluctuations is non-monotonic. This can also be 
seen in TablH For {K + K~)/{'i:~^ + vr"), the fluctuations decrease up to top SPS energy, 
afterwards very slowly increase with y/s. {p + P)/{tt^ + vr~) have negative dynamical 
fluctuations at SPS energies. At higher energies, their dynamical fluctuations might 
jump to positive values. The SPS energy might not high enough to lead to strongly 
out-of equilibrium phase transition. The situation is different at RHIC energies. On the 
other hand, the lattice QCD simulations show that the transition in the region of the 
phase diagram corresponding to the RHIC energy - and very early universe - is a cross 
over. This continuous phase transition is not strong enough to secure out-of-equilibrium. 
Furthermore, any anomalous phenomenon, like the dynamical fluctuations, should be 
washed out in the final state, as the transition is a smooth and continuous cross-over or 
second order. This might be the reason why although new state of matter should be 
created at RHIC energies |17) . there is no abrupt change in the dynamical fiuctuations 
with increasing y^. That the dynamical fluctuations smoothly increase with in agree- 
ment with the lattice simulations for Ug fluctuations [25; might support the conclusions 
that any anomalous phenomenon associated with non-equilibrium phase change likely 
would be washed out, if the phase change is smooth and does not cause out-of-equilibrium. 

In our analysis, we assume that the particle production is due to chemical equilib- 
rium processes controlling the final state, i.e., 7 = 1. That our models can very well 
reproduce the experimental measurements means that the equilibrium freeze out scenario 
is apparently correct, especially up to SPS energies. Nevertheless, energy scan down to 
^/s = 10 GeV turns to be a crucial step to verify the worthwhile behavior of particle 
production [14^ [T6] and now dynamical fluctuations. 

It is also worth extracting information about the role of different decay channels in the 
energy-dependence of dynamical fluctuations. It will be a further propose to study the 
effect of chemical non-equilibrium processes on event-by-event dynamical fluctuations. 
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